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I. INTRODUCTION 

Flows in networks are abundant: the flow of water 
in rivers and pipelines, the flow of current in electrical 
wires, the flow of passengers, the flow of cars through 
the network of roads [| are just a few examples. These 
flows can be characterized by the conservation of current 
- apart from sources and sinks in the network, the current 
is locally conserved. 

The maximal flows in capacitated networks are espe- 
cially interesting |2|. A capacitated network is a graph 
with a non-negative number called capacity c = c(e) as- 
signed to each edge e. Capacity measures the maximal 
flow that can pass through the edge. For each vertex, the 
current flowing in and out of it should be the same (Kir- 
choff's conservation law). This rule is only modified for 
the source vertices where the current enters the system, 
and the sink vertices through which the current leaves 
the network. Practically it means, that for two edges 
ei, e2 in series the combined capacity of the path (e\, e-i) 
is min(ci, C2); for two edges in parallel the combined ca- 
pacity is c\ + c<i. 

In this paper we mainly focus on rooted trees in which 
all leaf nodes are at the same distance from the root; 
such trees are often called perfect trees. We assume that 
the current enters through the root of the tree and it 
is discharged at the leaves (see Fig. [TJ. Generally for a 
path (ei,...,e n ), i.e., a set of edges in series, the ca- 
pacity of a path is the minimal capacity along the path: 
min(ci, . . . , c n ). This defines the flow between any two 
vertices in a tree. 

There is a huge engineering and mathematical liter- 
ature on flows in networks. One important goal is to 
find the maximal flow [3 03- C3- 0- E3| - although other is- 
sues have been also investigated, see e.g. fH \&L IqL fTot ITlT 
IT2I Il3l IbH ITU IHH . Here we concentrate on large net- 
works where the capacities are random variables chosen 
from the same distribution. We shall mostly investigate 
trees since the absence of circuits greatly simplifies the 
analysis. In Sec. [HI we investigate flows on deterministic 
rooted trees, namely on binary trees and then on b— ary 
trees. In Sec. IIIII we solve the maximal flow problem on 
hierarchical lattices which generalize rooted trees to in- 
clude loops. We then analyze random recursive trees in 
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FIG. 1: Illustration of a flow on a binary tree of three gener- 
ations. The width of the links represents their capacity, for 
instance c\ > C2. The current flows from top to bottom. The 
arrows represents the amount of current flowing in and out of 
the system. 



Sec. IIVI A brief discussion is given in Sec. 



II. DETERMINISTIC ROOTED TREES 

A rooted binary tree, also known as a Cayley tree of 
coordination number 3, has a root vertex in generation 
g = which is joined to 2 vertices in generation g = 1 
each of which is joined to 2 vertices in generation g = 2, 
so that there are a total 4 vertices in generation 2; gen- 
erally, the binary tree has 2 9 vertices in generation g 
(Fig. . Suppose that current can flow only in one di- 
rection, namely from generation g to g + 1. A capac- 
ity c = c(e) > is assigned randomly to each edge e. 
We shall assume that capacities are independently cho- 
sen from a distribution with density f(s) — —F'(s) where 



F(s) = Prob{c > s} . (1) 

For the binary tree with g generations, there are 2 9 
leaves (vertices in the last generation g). The capacity 
min(ci, . . . , c g ) of each path from the root to the leaf is 
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a random variable whose distribution is given by 

Prob{min(ci, . . . ,c g ) > s} = 
g 

JJ Prob(c fc >a) = [F(s)] 9 . 
k=i 

What is the probability distribution of the maximal flow 
among these 2 9 flows? What is the distribution of the 
total flux from the root to the leaves? These are the 
questions discussed below. 



A. Maximal flow 

Let M g = max(si, . . . , S2») be the maximal flow out of 
the 2 s flows to all leaves. Let 

P g (s)=Prob{M g <s} (2) 

be the cumulative distribution. Each path goes through 
one of the two edges issuing from the root (we denote 
their capacities by c\ and c-i as illustrated on Fig. 
Therefore the probability P g (s) is given by 

Prob{min(ci, M^) < s} x Probjmin^M^) < s } 
where M^\, M g \} 1 are the maximal flows in the corre- 



sponding daughter trees (whose roots are vertices from 
generation one). The daughter trees are statistically in- 
dependent and hence M^^, ±v± g _ x 
bution. Since 



r(1) , M {2 \ have the same distri- 



Prob{min(c,M g _i) > s} = F(s) [1 - P g -i(s)\ 
we arrive at a recurrence 

P g (s) = {l-F(s)+F(s)P g _ 1 (s)} 2 . 



(3) 



This recurrence could have been derived in many other 
ways. To have a smaller than s maximal flow in the tree, 
one needs, independently in both sub-trees, either the 
upmost capacity to be smaller than s {term 1 — F(s)}, 
or if it is larger than s, one needs the maximal flow in the 
daughter tree to be smaller than s {term F(s)P g -i(s)}. 
Iterating JSJ one gets 

Pi = {l-Ff 

Pi = [l-F + F{\-F) 2 } 2 

P 3 = |l-F + F[l-F + F(I-F) 2 ] 2 } 2 

etc. The expressions P g (s) become unwieldy but quickly 
approach a limiting distribution -Poo(s) = hirig->oo P g {s) 
which has a remarkably simple form 



P(s) 



1 - F(s) 



F(s) 
1 



for 1/2 < F(s) < 1 



for < F(s) < 1/2 



(4) 
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FIG. 2: Density p(s) of the maximal flow though a binary tree 
of g generations. The capacities were drawn independently 
from a flat distribution. Simulation results and analytic so- 
lution is depicted for finite g values. The g — > oo asymptotic 
limit is also shown. 



(Hereinafter we will omit the g — oo subscript.) The 
corresponding density function is 



p(s) 



{ 2f(s)[l-F(s)] 
[F( S )]3 





for F(s) > 1 /2 



for F(s) < 1/2 . 



(5) 



The most intriguing feature of this density function is 
the sharp cutoff at a*, which is given by F(s*) = 1/2. 
Note that the position of the cutoff s* does not depend 
on the details of the capacity distribution, only on the 
location it takes the value 1/2. In other words, the spe- 
cific values of the large capacities (s with F(s) < 1/2) do 
not matter for large trees. For instance, if the capacity 
distribution is bimodal, f(s) = (1 — a)S(s) + aS(b — s), 
then when a < 1/2 the maximal current is zero inde- 
pendent of b while for a > 1/2 the maximal current is b. 
This is essentially a branching process, where the process 
survives if the branching parameter is larger than 1/2. 

We now illustrate the behavior of the limiting distri- 
bution in two representative cases. As an example of a 
bounded capacity distribution, we choose the flat distri- 
bution with density f(s) = lif0<s<l and f(s) = 
when s > 1. Then 



p(s) 



) (T^F for s<1/2 

for s > 1/2 



(6) 



The sharp cutoff in the limiting distribution, and also the 
relatively slow convergence of the finite generation curves 
can be observed on Fig.|3 From @ one can compute any 
moment of the maximal flow. For instance, the average 
maximal flux is (s) — J ds sp(s) = 2 In 2 — 1. 
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As an example of an unbounded capacity distribution 
we take the exponential distribution, f(s) — e~ s . Then 



p(s) 



2e s (e s - 1) for s < In 2 



(7) 



for s > In 2 

and the average maximal flux is (s) = 1/2. 

B. Total flux 

Consider the total flux (more precisely the maximal 
possible total flux) from the root to the leaves. For the 
binary tree with one generation, the total flux from the 
root to the two leaves is <£>i = c\ + c 2 . Similarly for the 
binary tree with two generations 

$ 2 = min[ci, en + c i2 ] + min[c 2 , c 2 i + c 22 ] ■ 

Generally the total flux from the root to the bottom of 
the binary tree is found from 



$ g = min[ci, + min[c 2 , 



(8) 



where ® g li are the total fluxes in the correspond- 

ing daughter trees (whose roots are vertices from gener- 
ation one). 

Let Rg{q) = Prob($ s > q) and p g (q) = —R' g {q) be the 
corresponding density function of the total flux. As the 
daughter trees are statistically independent, equation (JHJ 
implies that 



Pg(Q) 



dx h g (x)hg(q — x) 



(9) 



where h g is the probability density of the total flux of a 
half tree 

h g (q) = -H' g (q) , H g (q) = Prob{min[c, $ 9 _i] > q} 

Since the minimum of two numbers are larger then q if 
and only if both numbers are larger than q, we have 



H g (q) = F(q) R g ^(q) 



(10) 



Setting Ro(q) = 0, Eqs. lf§|l- lfTn|> provide a recursive for- 
mula for p g {q) for any finite g. In the g — > oo limit we 
arrive at 



dx h(x)h(q — x) 



H(q) = F(q)R(q) . 



(11a) 
(lib) 



Equation Ijllajl shows that p(q) is a convolution. This 
suggests to employ the Laplace transform, which we de- 
note by a hat 5"(A) = J °° dqe~ Xq 3 r (g). We can recast 
()llajl into 



p(X) = [h(X)] 2 = [l-XH(X)f 



(12) 



and hence, with Ijllbfl . we have a closed equation for the 
Laplace transform of the flux density. 

First, let us discuss the properties of p(q) for a general 
capacity density f(q). Obviously, p(q) has a finite sup- 
port if and only if f(q) has a finite support. From l|lla|) 
and (|llb|l we find linear small q behavior p(q) — > qf 2 (0) 
in the generic case /(0) > 0. In Appendix |U we obtain 
the complete small q series of p(q). 

We now again consider two representative examples. 
For the exponential distribution of capacities, Eq. I|12|) 
becomes 



i + \p{\ + i) 

l + A 



(13) 



By iteration, one arrives at a formal solution 

>_/ 1 t A / 1 l + A / 1 \ 2 \ 2 \ 2 

P_ Vl + A + l + AV2 + A + 2 + AV3 + A + "'/ ) ) 

It seems hard to recast the Laplace transform into a 
more compact form, and to invert the Laplace transform. 
One can extract, however, the asymptotics of the flux 
density already from Eq. H13fl . The small q behavior of 
p(q) is encoded in the large A behavior of the Laplace 
transform. Using Eq. i|13|) one derives the large A expan- 
sion of the Laplace transform 



/5(A) = A- 



4A" 



2A" 



30A~ 



from which 



2 3 1 4 1 5 
P = q ~3 q + Y2 q + l q 



0. (14) 



This series can be also obtained from the general formula 
given in Appendix lAl 

If the large q behavior of p(q) were exponential, it 
would be encoded in the poles of its Laplace transform. 
We now show that the Laplace transform p(X) is an en- 
tire function, i.e., it has no poles or other singularities in 
the complex A plane. An apparent pole at A = —1 is not 
a pole — using the expansion 



P (X) = l-X(q)+0(X 2 ) 



(15) 



one finds that p{— 1) = (1 + (q)) 2 . Similarly one com- 
putes p(-2) = (1 + A{q) + 2(q 2 )) 2 and p(N) for other 
negative integers N and finds that all the apparent poles 
are regular points. The lack of poles means that p(x) 
decays faster than exponentially. 

To establish the decay law we use equation <|13|) to find 



/5(A) oc exp{«2~ A } 
for large negative A from which 



p(q) oc exp 



' In 2 



1 + In v 
In 2 



(16) 



(17) 
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FIG. 3: Simulation results for the density p(q) of the total 
flux though a binary tree of g generations. The capacities 
were drawn independently from an exponential distribution. 
The solid line represents the small q asymptotic solution given 
in Appendix 1X1 



as q — > oo. Thus p(x) decays faster than exponentially. 
The leading factorial asymptotic p oc e _?log 2 ? is inde- 
pendent on v, while the exponential correction depends 
on v whose determination requires the analysis of the full 
equation l|13fl . The asymptotic predictions (|14fl and i|17|) 
are in good agreement with simulation results (Fig. [31) ■ 
For the flat capacity distribution, Eq. Ijllb|) gives 



h(q) = (1 - q) p(q) + / dy p(y) for q < 1 (18) 

Jq 

and h{q) = for q > 1. Equation i|lla[l shows that p(q) 
vanishes outside the interval < q < 2, while inside the 
interval we utilize the fact that h(q) = when q > 1 and 
re-write (jllaf) as 



dx h(x)h(q — x) for < q < 1 



(19) 



dx h(x)h(q — x) for 1 < q < 2 



9-1 



The small q expansion of the flux density is found from 
Eqs. iP )l -iP )l to yield 



p(9J = q + q - ;9 -79 
6 6 



(20) 



in agreement with the general formula given in Appendix 
1X1 Similarly by expanding Eqs. dJ|-l© near the U P~ 
per cutoff we find that the flux density vanishes linearly 
according to 



p(q) -> ( 2 - q)p 2 > A* 



dx p(x) 



(21) 



The normalization requirement dx p(x) = 1 shows 
that p < 1, but to compute seems impossible with- 
out solving the entire problem. Numerically p ~ 0.4 



FIG. 4: Simulation results for the density p(q) of the total 
flux though a binary tree of g generations. The capacities 
were drawn independently from a flat distribution. The dis- 
tribution for the first generation is composed by two straight 
lines. The thick line represents the small q asymptotic solu- 
tion, valid for q < 1 and given in Appendix lAl 



and the asymptotic predictions (|2U)l and (|2T|) are in good 
agreement with simulation results (Fig.0J|. Interestingly, 
the total flux density converges remarkably fast in g, as 
opposed to the slow convergence of the maximal flow den- 
sity shown on Fig. [21 



C. 6-ary trees 

The above framework can be extended to other kinds 
of trees. For instance, for the &-ary tree (also known as a 
rooted Cayley tree with coordination number 6+1) each 
node (apart from the root) has one incoming link and b 
outgoing links. The analog of Eq. J2J reads 



P g (s) = {1 - F(s) + F(s)P g ^(s)} b 



(22) 



Therefore the stationary distribution P = P^ satisfies 
P = [1 - F + FP] b . The trivial root of this equation 
gives the limiting distribution for large s, 

P(s) = l when F(s)< 1-1/6 

while for F(s) > 1 — 1/6, the limiting distribution is given 
by a non-trivial root. For instance, for ternary trees 



P 



2F 2 -3F + V4:F - 3F 2 
W 2 



for F(s) > 2/3. In particular, for the ternary tree with a 
flat distribution of capacities 



P = 



1 + 3s 3s - 1 



(1- S )3 
1 



1 - s 



for s < 1/3 
for s > 1/3 
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FIG. 5: Illustration of a flow on a binary hierarchical lattices 
of one and two generations. The width of the links represents 
their capacity, and the arrows stand for the amount of current 
flowing in and out of the system. 



Further, the distribution of the total flux $ = $oo 
in the infinite rooted 6-ary tree is invariant under the 
transformation 

$ = min[ci, $] + . . . + min[cb, $] . (23) 

Therefore the flux density is a (generalized) convolution 

ri ri b 

p(q) = ... Y[dxjh(xj) 8{q-xx- ...-x b ) (24) 
Jo Jo - =1 

and the Laplace transform is /3(A) = [ft(A)] b . As an ex- 
ample, let us consider the exponential distribution. Then 
the Laplace transform of the flux density becomes 



/5(A) 



l + Xp(X + l) 



i b 



1 + A 



The small q expansion is 



P 
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6-1 



0(q b ) 



(25) 



(26) 



(6-1)! 

while when q — > oo the flux density decays according to 
p(q) oc e - qlo ^ q (27) 

III. HIERARCHICAL LATTICES 

Hierarchical lattices represent a simple generalization 
of rooted trees. These lattices have loops but are still 
tractable (see e.g. [HE El). They mimic finite di- 
mensional lattices. 

A 6-ary hierarchical lattice of g generations is com- 
posed by two 6-ary trees of g generations merged at the 
leaves. The current enters the system at one root, and 
leaves it at the other. In Fig|S]we show the binary hier- 
archical lattices of one and two generations. 

Denoting again the flux distribution of a half of the 
lattice by H g we get 



In contrast to Eq. I|10(l we have F 2 since the capacities 
of two edges attached to the daughter lattice with g — 1 
generations should be smaller than q. 

The total flux is given by a 6-fold convolution analo- 
gous to (|24|l . which again simplifies after performing the 
Laplace transform 



p g (\) = {h g (\)] b = [l-\H g (\)] b 



(29) 



A finite g distribution can be obtained by starting with 
Ro(q) = and iterating the above equations. 

The limiting flux distribution is determined by solving 



H(q) = [F(q)fR(q) 
p(\) = [l-\H(\)] b 
In the generic case of /(0) > we get 
[2/(0)]" 



P(8) 



0(q b 



(30a) 
(30b) 

(31) 



(6-i)r 

in the small q limit. 

The large q behavior is less universal. Let us consider 
again two representative cases exemplifying unbounded 
and bounded capacity distributions. For the exponential 
capacity distribution, Eqs. H30a[) - (|30b(l give 

+ xp(\ ■ > " 



P(A) 



2 + A 



(32) 



Using this expression we extract the asymptotic behav- 
iors. We find that the small q expansion agrees with 
the general asymptotic l|31l) while when q — > oo the flux 
density decays according to 

p{q) cx e - 2qlos » q 

These predictions agree with simulation results (Fig.[f)J). 

For the flat capacity distribution we arrive [by differ- 
entiating Eq. I|30al) ] at an integral equation 

r-2 



h(q) = (1 - q) 2 p{q) + 2(1 - q) / dy p(y) 



(33) 



for q < 1. Another relation between functions p(q) and 
h(q) is (in the case of the binary hierarchical lattice) given 
by Eq. (fTU]b ■ 

To extract the behavior of p{q) in the q — > 2 limit we 
write q = 2 — e, x = 1 — ey and re- write Eq. I|19l) as 



P(?) 



: f dyh[l-ey]h[l-e(l-y)} 
Jo 



(34) 



From Eq. we get 

h[l-ez] = 2epz + 0(e 2 ) (35) 

where again p — J 2 dx p(x). Plugging (f^3f) into Eq. 
and computing the integral we find that p(q) vanishes 
cubically near the upper cutoff, 



p(q) = -p 2 (2-qf + 0[(2-q) i ] 



(36) 



H g (q) = [F(q)} 2 R g _ 1 (q) 



(28) 



Generally for the 6-ary hierarchical lattice p oc (2 — q) b+1 
near the upper cutoff. 
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which is also solvable. Asymptotically 



FIG. 6: Total flux density p(q) for binary hierarchical lattice 
in the infinite size limit (g — * oo). The capacities were drawn 
from a flat and from an exponential distribution. The curves 
were obtained by iterating 1281 , and recovered by the small q 
asymptotic series given in Appendix 1X1 



IV. RANDOM RECURSIVE TREES 

In the previous sections we have studied flows on reg- 
ular graphs where all node degrees (apart from the root 
node) were the same. It is harder to deal with arbitrary 
random graphs but at least one class of random trees, 
the so-called random recursive trees, are tractable. First 
we overview the general properties of these random trees, 
and then we shall discuss flows on these trees. 

A random recursive tree is generated as follows: One 
starts with a root and adds nodes one by one so that each 
newly-introduced node is linked to a randomly selected 
existing node. Random recursive trees which have been 
studied in great detail (see e.g. the survey [2(| and more 
recent articles HH H S || ) . 

Let N be the total number of nodes and Dj (A) be the 
number of nodes on distance j from the root. By defi- 
nition, Do(N) = 1, since the root is on distance from 
itself. Quantities Dj(N) with 1 < j < A are random, 
and e.g. the number of nodes on distance one from the 
root evolves according to 



Di(N) + l prob A" 1 
Di(iV) prob 1-N- 1 



(37) 



Thus the average Di(iV) = (D\ (N)) satisfies 

V,(N+ 1) = D, (AT) + 1 

whose solution is Di(N) = J2i<j<N-i J -1, Similarly us- 
ing (|37|) one establishes a simple recurrence for the vari- 
ance Vi(A) = (Df(AO) - (£>i(A)) 2 , viz. 



Vi(JV + l) = Vi(N) + 



N 



1 



2) 1 (A)=lnA + 7 + 0( — 



V 1 (A)=lnA + 7 -^ + 0(l 



(38a) 
(38b) 



where 7 = 0.5772 ... is Euler's constant. 

Thus fluctuations are asymptotically negligible [2{| 
and we can focus on the averages. One easily establishes 
the exact recurrence for the averages 

■Dj (N + 1) = ^ (AT) + ^D J --i (AO • (39) 

In the N — > 00 limit, this recurrence reduces to a dif- 
ferential equation 4^ = A -1 Dj-i whose solution reads 



(In Ny 



(40) 



Thus the root is linked to approximately In A nodes 
each of which in turn is linked on average to i In A nodes. 
What is important is that each node in the first gener- 
ation (on distance one from the root) is linked to many 
nodes in the second generation; only a few nodes in the 
first generation have a finite (not diverging with A) num- 
ber of links to the nodes of the second generation. This 
observation allows us to compute the maximal flow in the 
leading order. Indeed, asymptotically it is the maximal 
capacity among all D\(N) capacities from the root to the 
first generation nodes. 

The above argument shows that for the flat distribu- 
tion of capacities the density distribution for the maxi- 
mal flow approaches a stationary limit Poo{s) = 5(s — 1). 
Essentially the same behavior occurs for any bounded 
capacity distribution: Poo(s) = S(s — s max ). For un- 
bounded capacity distributions, the density distribution 
for the maximal flow does not approach a stationary 
limit. For instance, for the exponential distribution we 
approximately have 



-s\lniV 



PN (s) =lnAe- s (l-e- s ) 
implying that on average the maximal flow is 



{s)n 



In AT 



(41) 



(42) 



The computation of the total flux is even simpler. For 
an arbitrary capacity distribution, the total flux from 
root to leaves is 



N 



(c) I11N 



(43) 



Fluctuations around the average are theoretically negli- 
gible in the A — > 00 limit, but since they scale as Vhi A, 
see Eq. (|38b|l , fluctuations become negligible only for im- 
mensly large A. 
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V. DISCUSSION 

Trees represent the least interconnected graphs while 
complete graphs are the most connected. Flows in com- 
plete graphs exhibit simple asymptotic behavior. The 
maximal flux from an arbitrarily chosen source node to 
an arbitrarily chosen sink node scales as 



(s) N = (c)N 



(44) 



and fluctuations about this average are asymptotically 
negligible. To show this we first notice that for a com- 
plete graph with N + 1 nodes, the flow from node 
(source) to node N (sink) cannot exceed J2i<j<N ^ c oj 
which is asymptotically (c)N (fluctuations are of the or- 
der of \/N). This is the upper bound and from a par- 
ticular node j it will be impossible to transmit flow CQj 
directly to the sink if CQj > Cjjy However the node j is 
connected to nodes 1 , . . . , j — l,j + 1 , . . . N — 1 and in 
the large N limit it will be possible to find the way to 
transfer the flow via those side routs. 

For regular lattices, the problem seems analytically 
intractable. The simplest set-up is a two-dimensional 
square lattice with additional restriction that the flow is 
biased along the diagonal. More precisely, the root is 
(0,0) and from every site with i,j > the flow is 
to (i + 1, j) and + 1). The capacities are again cho- 
sen independently from the same distribution. A similar 
problem of directed polymers on the same lattice with 
random 'energies' assigned to nodes has been solved by 
Krug and Halpin-Healy [26|. The difficulty in our case 
arises from the splitting of current at nodes. 

A related problem is the minimal cost flows on net- 
works ■ In these models there is a cost associated with 
each edge and the cost of a flow is the sum of the costs of 
the edges it flows through. This problem was considered 
with random costs drawn from a common distribution 
0, 0| . When the costs r are taken from an exponential 
distribution exp(ar), the a — > oo limit is called the strong 
disorder limit. In this case the cost of a flow is dominated 
by its single most expensive edge, and the search for the 
cheapest flow becomes equivalent to our present study on 
trees: the search for the maximal flow. 



where • W denotes the i-th derivative of a function. The 
integrated probabilities then can be written as 



^ /(*-P (Q) , -pO-i)(0) 

F (q) = -2^ — j\ — 9 , R{q) = -2_^ — j\ — i 



i=0 



3=0 



where formally we denote /( ^(0) = = — 1- 

Then from ljllb|) we have 



n=0 i=0 



/(»-!) (Q)p("-<-l) (Q) 

i\(n — %)\ 



from which we obtain 



-&W(0) 



i=0 



i!(n+ 1 - i)\ 



On the other hand, from Ijllajl 



p( fe )(o) = ^/i( J '- 1 ) (0)^-^(0) 



(A3) 



Hence, with (|A2|I and (|A3() . we explicitly expressed 
(0) with lower derivatives of p, which recursion can be 
easily implemented on a computer. The limiting distri- 
bution is compared to simulation results in Fig. [3] and 0] 
(for the flat capacity distribution the result is valid only 
for q < 1, as the radius of convergence is one for the flat 
distribution itself). It is also straightforward to obtain 
the leading terms of the asymptotic series by hand 



P(q) = qf(0) + 9 2 /(0)[/ 2 (0) + f'(0)} + 0(q 3 
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APPENDIX A: SMALL q SERIES OF p(q) 

Assume that we can write f(q) and p(q) as a series 
around q = 

rf.Hf^ (A1) 

i=0 ' j=0 J ' 



For binary hierarchical lattices the calculation is simi- 
lar. Equation i|A3() is still valid, while one should replace 
llAHl by 

^ B .(, +lfE a50sa!!, 

i,j,l>0 

where i + j + l = n + lis required in the sum. 

The above calculations can be straightforwardly ex- 
tended to &-ary trees and 6-ary hierarchical lattices. 
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